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ABSTRACT 

Millimeter-wavelength VLBI observations of the supermassive black holes in Sgr A* 
and M87 by the Event Horizon Telescope could potentially trace the dynamics of 
ejected plasma blobs in real time. We demonstrate that the trajectory and tidal 
stretching of these blobs can be used to test general relativity and set new constraints 
on the mass and spin of these black holes. 


1 INTRODUCTION 

The planned Event Horizon Telescope (EHT(Qwill possess 
angular resolution comparable to the Schwarzschild radius 
of the supermassive black holes (SMBHs), Sgr A* and the 
one at the center of M87, and temporal resolution on min¬ 


utes timescales (Johnson et al. 20141. This is expected to 


open a new avenue for studying a multitude of transient 
phenomenae under extreme gravity. 

Sgr A* is known to exhibit variability with tens of min¬ 
utes timescale corresponding to accretion disk activity at the 
innermost stable circular orbit (ISCO) (Genzel et al.|[2003 


[Johnson et al.|[2014[ ). Here we study a hypothetical class of 
short timescale events corresponding to plasma blobs ejected 
near the ISCO radius. Although such blobs were never ob¬ 
served from a supermassive black hole, they may exist based 
on the analogy with microquasars, which are known to pro¬ 
pel blobs at relativistic speeds (Mirabel et al. 1992 Mirabel 
[200^|2004| ). 

The second target of the EHT is the supermassive black 
hole at the center of the elliptical galaxy M87. In contrast 
to Sgr A*, M87 possesses a jet, and it is likely that blobs 
are ejected along the jet’s symmetry axis. 

In this Letter, we demonstrate that if ejected plasma 
blobs were detected, one could use their dynamics to probe 
the spacetime around the black holes. Furthermore, if the 
mass and spin of a given black hole are known, one can use 
observations of the blob’s dynamics to test general relativ¬ 
ity or infer the presence of non-gravitational sources such as 
gas pressure or magnetic stress. These constraints would be 


complimentary to constraints from pulsar timing | Pfahl and 


Loeb||2004| [Cordes et al.||2002| [Kramer et al.||2004[ jPsaltis 

and Johannsen|2011[ Liu et al.|2012 1 or observations of the 


black hole shadow |[Lu et al. 2014 
hannsen|2012 1. 


Psaltis et al. 2014 Jo- 


There are two elements of dynamical information: the 
trajectory of the blob’s center of mass, and its lateral expan- 


^ http://www.eventhorizontelescope.org/ 


sion. Both can be used to independently constraint the black 
hole’s spacetime. We discuss the former in §2 and §3, and 
the later in §4. Throughout the discussion, we will assume 
general relativity. Deviations from our results would indi¬ 
cate the presence of non-gravitational forces or corrections 
to the theory of gravity. We use units where G = c = 1, and 
the conversion from these units to physical units is given in 
Tabled! 


2 CENTER OF MASS MOTION 


First we consider the motion of the blob’s center of mass 
(COM). If the blob is ejected above the escape speed from 
the ISCO radius, Risco, its azimuthal velocity will be negli¬ 
gibly small at r ^ Risco, so we focus our discussion on the 
radial equation of motion. For a Schwarzschild black hole 
I Chandrasekhar[|1983 1, 
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where M is the black hole mass, e the energy per unit rest 
mass of the blob, r the black hole-blob distance, t the coordi¬ 
nate time, and r the blob’s proper time. These two equations 
can be solved for dt/dr and integrated to obtain the coor¬ 
dinate time as a function of the orbital radius of the blob’s 
COM, 
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If the blob is ejected out of a Kerr black hole, a similar set 
of equations can be solved to obtain its COM motion in the 
equatorial plane, 
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where a is the black hole’s spin parameter and A(r) = — 

2Mr + a^. In general, there is no reason for the blob to be 
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Black hole mass 

Distance 

Time 

Space 

Angle 

Sgr A* 

(4.31 ± 0.36) X 10® Mq 

7.94 ± 0.42 kpc 

21 s 

0.043 AU 

5.3 pas 

M87 

(6.6 ±0.4) X 10®M0 

16.7 ± 0.9 Mpc 

9 hr 

65 AU 

3.9 pas 


Table 1. The conversion of black hole mass, M, to units of time, space, and angular size on the sky for Sgr A* and M87 ( [Gillessen et| 
|al.|2009l |Ghez et al.|200^ [Walsh et al.|20T^ [Eisenhauer et al.|2003| [Gebhardt et al.|2011| | , for G = c = 1. 



Figure 1. The radial motion of blobs with e = 2 in the equitorial 
plane of a black hole with a = 0 (solid line) and a = 0.999 (dotted 
line). 


ejected in the equatorial plane of the black hole, and in fact 
blobs should preferentially be ejected along the spin axis. 
But, as shown in Figure the effect of the black hole spin 
is weak. At t = lOM, the trajectory of a blob with e = 2 
launched from an a = 0.999 black hole is only 0.36M apart 
from one launched from an a = 0 black hole. 


3 SIMULATED OBSERVATIONS 


In simulating what would be seen by radio interferometers, 
we project the COM motion of the blob to the sky plane 


far from the black hole. We utilize the geokerr code (Dex¬ 


ter and Agol 20091 to trace rays from the observer plane 


located at inhnity to the position of the blob. The coordi¬ 
nates {x, y) parameterize positions in this observer plane. 
The Fourier transform of this plane yields the visibility of a 
radio interferometer. 

The blob itself is modeled as a small sphere that is emit¬ 
ting isotropically in its rest frame. The result for blobs with 
velocity vectors at angles 6 = Q and 6 = n/?> away from the 
observer are presented in Figure For a blob moving along 
the 6 = 0 axis, the image is briefly lensed into a ring with ra¬ 
dius Rrino ~ 5M. Previous calculations by|Johannsen and| 


Psaltis (20101 showed that the eccentricity of this ring is 
not sensitive to the spin of the black hole (except for a « 1), 
but is very sensitive to the black hole’s quadrupole moment. 
Thus, if detected, the ring can be used as a test of the no¬ 
hair theorem. As the ring only appears when the blob is 
still close to the black hole, its lifetime is short (~ 40M for 
a blob with e = 10, but longer for slower moving blobs). It is 
therefore necessary to have temporal resolutions on minutes 
timescale to detect the ring. 

In addition, if the motion is fast enough and is launched 
at a small angle relative to the observer, the apparent tra¬ 
jectory can appear superluminal (e.g. Rees||19^|. Close to 
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Figure 2. Blobs with e = 10 and radius M launched with 6 = 0 
and 6 = 7r/8 as seen in the observer plane with the black hole 
located at (0,0). The observer’s time axis (in units of M) is indi¬ 
cated by the color bar. For a blob moving with 6 = 0, the image 
is briefly lensed into a ring. The eccentricity of this ring can be 
used to test the no-hair theorem. 


the black hole, this apparent superluminal motion will be 
obscured by the bright photon ring. Thus, the detection of 
superluminal motion will require either waiting for the ring 
to dim or a manual removal of the ring. 

The projected distance as a function of observed times, 
shown in Figure can be compared with observations to 
determine the presence of non-gravitational forces (e.g. due 
to magnetic fields or hydrodynamic friction on background 
gas). In addition, it can be used to constrain gravitational 
theories that predict changes on the orbit of test particles 
close to a black hole (e.g. Giddings|2014 1 . 


4 TIDAL EFFECTS 

If the forces holding the blob together are much smaller 
than the tidal gravitational forces, the blob will be tidally 
sheared. The magnitude of this tidal shear depends on the 
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Figure 3. The projected position of blobs with e = 10 launched 
at a variety of angles versus observer time. 



Figure 4. The growth factor of the blob radius due to gravita¬ 
tional tide as a function of distance from the black hole for a blob 
moving with negligible angular momentum. The blob’s specific 
energy is e = 1.0001, 1.001, 1.01, and 10 for the solid, dashed, 
dotted lines, and dot-dashed lines, respectively. 


black hole’s mass and spin, and thus can be used to probe the 
black hole metric. Under the approximation that the force 
per unit mass keeping the blob together is <C {2MR/P), 
where R is the radius of the blob, the elements of the blob 
can be treated as if they are moving along geodesics. 

If the blob is small, we can define the geodesic deviation 
vector between the geodesic followed by the particle at 
the center of the blob and the different geodesic followed by 
particles at the blob’s edge by, 




dx°‘ 

ds 


(4) 


where s is the parameter indexing neighboring geodesics. We 
can calculate the rate of change of with respect to the 
affine parameter of the geodesic, 


= uNpC - 

= . 
where we have used the identity 
uf^VpC = 


Poisson|[2004 1, 

Vpu°‘ 


(5) 

( 6 ) 

(7) 


which is valid for geodesic deviation vectors. Writing explic¬ 
itly, 

, ( 8 ) 

yields 


d pdu'^ 

dr dx^ 


(9) 


The four velocity of a blob ejected from a Schwarzschild 
black hole with negligible angular momentum is: 


u 


a 




(l-e2),0,0 


( 10 ) 


For relative motion between particles at the center of the 
blob and particles at the edge of the blob in the radial di¬ 
rection: 


r = (o,i?,o,o)“. 


Plugging equation (111 into equation gives: 


1 dR 

rIx 
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( 12 ) 


Note that substituting t for A in equation ( |12[ ), then taking 
a derivative with respect to t with M/r —^ oo reproduces the 
tidal acceleration of Newtonian gravity: audai ~ MR/P. 
Substituting the orbital radius r in place of A in equa¬ 


tion (121 and integrating, we get: 


dR' _ r 
Jro J ro 


Mdr' 


r'2 (_l + e2 + ^) 


(13) 


where Ro r is the initial size of the blob and ro the 
starting orbital radius of the blob. Assuming that the blob 
is ejected from the ISCO radius, ro = 6M for a = 0, we 
obtain: 


_R _ r (-2-b3e^)r 
Ro 6M -I- 3(e2 — l)r 


(14) 


This change in radius is in principle observable, and can 
therefore be used to find the mass of the black hole if e 
is inferred from the COM trajectory. The constant e can 
be inferred far away from the black hole where it obeys 
e = l/-\/l — WooAf, where vcom is the COM velocity of the 
blob Sbt r M. Figure]^ shows the radial growth factor for 
blobs with specific energy e = 1.0001, 1.001, 1.01, and 10. 
Because blobs of smaller e spend more time close to the black 
hole, the tidal effect is larger the closer e is to unity. In the 
case of e ~ 1, one can get a growth factor of R/Ro ~ 10 
at r = lOOOM. This is a change that is observable by the 
EHT. In general, one can also compute the relative motion 
between the center and the edge of the blob in the 4> and 6 
direction via an analogous calculation. 

We can extend this calculation to the case of a spinning 
black hole with a blob moving radially in the equatorial 
plane. For this configuration, the relevant components of w“ 
are, 

(15) 


(16) 

Again we adopt. 


u = ^ I ^ + o + 


2pM 


I 2^ 2M a2g2 A 

= * 'e2 H-5-(ae)2 -|---- . 


r = (0,i?,0,0)“ . (17) 


Performing an analogous calculation as in the a = 0 case. 
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Figure 5. The growth factor of the blob radius as a function 
of distance from a spinning black hole for a blob trajectory with 
a negligible angular momentum. The black hole’s spin is a = 
0, 0.5, and 1 for the solid, dashed, and dotted lines, respectively. 
The blob energy is e = 1.2 for all curves. 

we obtain, 

_R _ 

Ro 

Rosja" (-3 + 4e2) M + 36 (-2 + Se^) 

[2M + (—1 + e^) r] + Sa^ [—r + e^(2M + r)] 

(18) 

If the mass of the black hole and the blob energy e are 
known, this equation can be used to measure the spin of 
the black hole. Figure shows the growth factor RjRo for 
blobs with dimensionless spin parameter a = 0, 0.5, and 1. 
The effect of spin is weak, and its measurement would be 
challenging. 


Cordes, J. M. et al. 2002, New AR, 48, 1413 
Dexter J. and Agol, E. 2009, apj, 696, 1616 
Eisenhauer F. et al. 2003, ApJ, 597, L121 
Gebhardt K. et al. 2011, ApJ, 729, 119G 
Genzel, R. et al. 2003, Nature, 425, 934 
Ghez, A. M. et al. 2008, ApJ, 689, 1044 
Giddings, S. B. 2014, PhysRevD, 90, 12 
Gillessen, S. et al. 2009, ApJ, 692, 1075 
Johannsen T. 2012, RASP, 124, 1133 
Johannsen T. and Psaltis D. 2010, ApJ, 718, 446J 
Johnson M. D. et al. 2014, ApJ 
Kramer, M. et al. 2004, New AR, 48, 993 
Liu K. et al. 2012, ApJ, 747, 1 
Lu R. S. et al. 2014, ApJ, 788, 120 
Pfahl E., Loeb A. 2004, ApJ, 615, 253 
Poisson, E. 2004, A Relativists’s Toolkit, Gambridge Uni¬ 
versity Press, Gambridge 
Mirabel, I. F. et al. 1992, Nature, 358, 215M 
Mirabel, I. F. 2002, lAU Symposium, 214, 201 
Mirabel, 1. F. 2004, ESA Special Publication, 552, 175 
Psaltis D. et al. 2014, arxiv: 1411.454 
Psaltis D., Johannsen T. 2011, JPhGS, 283, 2030P 
Rees M. J. 1966, Nature, 211, 468R 
Walsh, J. L. et al. 2013, ApJ, 770, 86W 


5 CONCLUSION 


We have shown that observations of ejected plasma blobs 
from the supermassive black holes Sgr A* and M87, can 
be used to constrain the spacetime near these black holes. 
There are two pieces of information that can be obtained 
from these observations: the blob’s trajectory and the tidal 
effects on the blob’s shape. 

The trajectory of the blob can be used to limit the pres¬ 
ence of non-gravitational forces around the black hole or to 
constrain theories of gravity that predict anomalies in the 
orbit of test particles in the vicinity of black holes (e.g. Gid 
dings 20141. If a photon ring is detected, its eccentricity 
could be used as a test of the no-hair theorem. Eurther- 
more, observations of the tidal stretching of the ejected blob 
can be used to determine both the mass and spin parameter 
of the black hole. 
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